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1 Introduction and definitions 



1.1 Introduction 

This paper is the continuation of our study of surface phase transitions of 
the semi-infinite Potts model m (to be referred as paper I). 

Semi-infinite models exhibit a great variety of critical phenomena and we 
refer the reader to Ref. for a review on this subject. 

We consider, the g-states Potts model on the half-infinite lattice with 
bulk coupling constant J and surface coupling constant K (see below 
for the definition of the Hamiltonian) . 

Besides its popularity, this model presents very interesting features. Na- 
mely, in the many component limit g — > oo, the mean field theory yields by 
looking at the behavior of a bulk and a surface order parameter, and after a 
suitable rescaling i.e. by taking the inverse temperature (3 = Ing, the phase 
diagram shown in Figure 1 [Tfj . 




FIGURE 1: Mean field diagram borrowed from Ref. |27| . 

In region (I) (respectively (IV)) the bulk spins and the surface spins are 
disordered (respectively ordered). In region (II) the surface spins are ordered 
while the bulk spins are disordered. The region (III) called new low temper- 
ature phase [2Zj corresponds to disordered surface spins and ordered bulk 
spins: this phase which is also predicted by renormalization group scheme, 
actually does not appear in the Ising case [T7j. On the separating line between 
(I) and (IV) an ordinary transition occurs whereas so-called extraordinary 
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phase transitions take place on the separating hnes (I)-(III) and (II)-(IV). 
Finally, on the two remaining separation lines (I)- (II) and (III)- (IV), surface 
phase transitions arise. 

In paper I, we studied the high bulk temperature regime showing that the 
first surface phase transition between a disordered and an ordered surface 
while the bulk is disordered holds whenever e^*^ — 1 < q^^'^ and q is large 
enough. 

We are here concerned with the more interesting situation in which the 
bulk is ordered. We prove that the second surface transition between the 
new low temperature phase and the ordered phase actually occurs whenever 
qPJ — 1 > gi/*^^ again for large values of q. 

The results are based on the analysis of the induced effect of the bulk on 
the surface. Intuitively, this effect might be viewed as an external magnetic 
field. When the bulk is completely ordered (a situation that can be obtained 
by letting the coupling constant between bulk sites tends to infinity) the 
system reduces to Potts model in dimension d — 1 with coupling constant 
K submitted to a magnetic field of strength J . Such a model is known to 
undergo a order-disordered phase transition near the line (3J{d — 1) -|- (3K = 
In g ||5J. We control here this effect up to e^-^ — 1 > q^^'^ by a suitable study of 
a surface free energy and its derivative with respect to the surface coupling 
constant, which contains the thermodynamic of the surface phase transition 
under consideration. 

The technical tools involved in the analysis are the Fortuin-Kasteleyn 
representation (TH], cluster-expansion [121 1201 CDl I2H1 , Pirogov-Sinai theory 
as already in paper I, but in addition Alexander's duality [H 1211 1211 EOj ■ 

The use of Fortuin-Kasteleyn representation is two-fold. It provides a 
uniform formulation of Ising/Potts/percolation models for which much (but 
not all) of the physical theory are best implemented (see ^3] for a recent 
review). It can be defined for a wide class of model, making results easier to 
extend (see e.g. [221EnilHl)- This representation appears in Subsection 2.1 
and at the beginning of Subsection 2.2 to express both partition functions 
{Z and Q) entering in the definition of the surface free energy in terms of 
random cluster model. 

Alexander's duality is a transformation that associates to a subcomplex 
X of a cell-complex K the Poincare dual complex []K\X]* of its complement. 
Alexander's Theorem provides dualities relations between the cells numbers 
and Betti numbers of X and those of []K\X]* (see e.g. [Tll2ni)- FK measures 
on lattices are usually expressed in terms of the above numbers for a suitably 



3 



chosen cell-complex associated to the lattice under consideration. Alexan- 
der's duality provides thus a transformation on FK configurations (and FK 
measures) [2j. In the case of the Ising/Potts models this transformation is 
in fact the counterpart of the Krammers-Wannier duality (or its generaliza- 
tions P^l^l^ ): applying it after FK gives the same result than using first 
Krammers-Wannier duality and then taking FK representation [201 E] ■ We 
use Alexander's duality first in Subsection 2.2. It allows to write the bulk 
partition function (Q) as a system of a gas of polymers interacting through 
hard-core exclusion potential. The important fact is that the activities of 
polymers can be controlled for the values of parameters under consideration. 
This partition function can then be exponentiated by standard cluster ex- 
pansion. This duality appears again in Subsection 2.3 to obtain a suitable 
expression of the partition functions {Z ). 

Cluster expansion is used again in Subsection 2.3 to express the ratio 
Z/Q as a partition function of a system called Hydra model (different from 
that of paper I) invariant under horizontal translations. 

Pirogov-Sinai theory, the well-known theory developed for translation in- 
variant systems, is then implemented in Section 3 for the study of this system. 
Again cluster expansion enters in the game and the needed Peierls condition 
is proven in Appendix. 

The above description gives the organization of the paper. We end this 
introduction with the main definitions and a statement about the surface 
phase transition. 

1.2 Definitions 

Consider a ferromagnetic Potts model on the semi-infinite lattice L = Z'^"^ x 
Z+ of dimension d > 3. At each site i = G L, with ia & for 

a = 1, ...,d—l and id G Z+, there is a spin variable cXi taking its values in the 
set Q = {0, 1, . . . , q — l}. We let d{i,j) = max^^i^ ..^^ \ia — ja\ be the distance 
between two sites, d{i,Q) = minjizn d{i, j) be the distance between the site 
i and a subset Q C h, and d{fl,fl') = minifzn j d{i, j) be the distance 
between two subsets of L. The Hamiltonian of the system is given by 



where the sum runs over nearest neighbor pairs (i.e. at Euclidean dis- 
tance d-E{i,j) = 1) of a finite subset C L, and 6 is the Kronecker sym- 




(1.1) 
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bol: (5((7j,(7j) = 1 if (Tj = aj, and otherwise. The couphng constants K^j 
can take two values according both i and j belong to the boundary layer 
Lo = {i e L I id = 0}, or not: 

A«=l 'Pn (1-2) 
^ [ J > otherwise ^ ^ 

The partition function is defined by: 

Z-(0)^5]e-'^Vn (1-3) 

Here the sum is over configurations ctq G Qp, (5 is the inverse temperature, 
and Xq is a characteristic function giving the boundary conditions. In par- 
ticular, we will consider the following boundary conditions: 

• the ordered boundary condition: Xn — Wi^an H^i^ 0); where the bound- 
ary of fl is the set of sites of Q at distance one to its complement 
an = {i e n : ci(i,L\0) = 1}. 

• the ordered boundary condition in the bulk and free boundary condition 
on the surface: Xn — Iliedbn^i'^i^ 0)' where dbQ — dQ f] (L \ Lq). 

Let us now consider the finite box 

f2 = {i e L I max \ia\ < L, ; < < M} 

a=l,...,d—l 

its projection E = fl Lq = {i G ^2 | = 0} on the boundary layer and its 
bulk part A = Q\E ^ {i e Q \ 1 < id < M}. 
The ordered surface free energy, is defined by 

1 Z°(fl) 
Qo = — hm ■— — lim In — — (1.4) 

^° L-.0O |E| M-»oo Q°(A) ^ ^ 

Here |S| = (2L + 1)''^^ is the number of lattice sites in E, and (5°(A) is the 
following bulk partition function: 

g°(A) = ^exp{/?J 5(<7„<7,)} n^i^O) 

{i,j)cA ieOA 

where the sum is over configurations ua G Q^. The surface free energy does 
not depend on the boundary condition on the surface, in particular one can 
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replace Z°{fl) by Z°^{Q) in (jl.4|) . The partial derivative of the surface free 
energy with respect to I3K represents the mean surface energy. As a result 
of this paper we get for q large and q^^^ < e^'^ — 1 < g that the mean surface 

energy g^go is discontinuous near PK = lnil+( -^pj-^j j j . 

Namely, let ( ■ y denote the infinite volume expectation corresponding to 
the boundary condition p: 

{fy{pj,PK)= lim y /e-^^Xc 

defined for local observable / and let be defined by ()3.8|) below. As 
a consequence of our main result (Theorem 13.51 in Section 3), we have the 
following 

Corollary 1.1 Assume that q^^^ < e^'^ — 1 < q and q is large enough, then 
there exists a unique value Kt{P, J, q, d) such that for any n.n. pair ij of the 
surface or between the surface and the first layer 

{6{ai,aj)yf{(3J,(3K)<Oie-^) for K < 
{5{a.,,a,)Y{(3J,l3K)>l-0{e-^) for K > 

In that theorem the ratios of the partition functions entering in the defini- 
tion of the surface free energy g^ (with both Z°{VL) and Z°^(f2)) are expressed 
in terms of partition functions of gas of polymers interacting through a two- 
body hard-core exclusion potential. For g^/*^ < e^'^ — 1 < q and q large, the 
associated activities are small according the values of K namely for K > Kt 
with the ordered boundary condition and for K < Kf with the ordered-free 
boundary condition. The system is then controlled by convergent cluster 
expansion. 



2 Random cluster models and Hydra model 

2.1 The Fortuin— Kasteleyn (FK) representation 

By using the expansion g/^-f^ij^C'^^.'^j) = 1 + (e^^»j — l)5{ai,aj), we obtain the 
Fortuin-Kasteleyn representation [16 of the partition function: 

XcB{n) (i,j)€X 



6 



where B{Q,) = {{i,j) : i E fl, j & fl} is the set of bonds with both endpoints 
belonging to and N^(X) is the number of connected components (regard- 
ing an isolated site i G ^2 as a component) of a given X C B{Q). These 
numbers depend on the considered boundary condition; introducing S{X) as 
the set of sites that belong to some bond of X and C{X\V) as the number 
of connected components (single sites are not included) of X that do not 
intersect the set of sites V, they are given by: 

iV°(X) = \Q\ - \S{X) U dQ\ + C{X\dQ) 
K\X) = \n\- \S{X) UdbQ\ + C{X\dbQ) 

Hereafter \E\ denotes the number of elements of the set E. 
We introduce the parameters 

ln(e^-^ - 1) 
~ hi 

_ In(e^'^-l) ^^-^^ 

Pb = ] 

m q 

and let X^ ^ X D S(Lo), Xb ^ X \ Xs, to get 

ZP{n)= ^ qm+Pb\x,\+Nf,ix) ^2.3) 

X<ZB{Q) 

The ground state diagram of this system is analogous to the diagram of 
Figure 1, by replacing J by Pb and K by (3s (see paper I). 

For the bulk partition function Q°{A), one find that the FK representation 
reads 

g°(A) = ^ qMY\+Nl{Y) ^ ^mA)] ^ ^-mA)\Y\+Nl{Y) ^2A) 
Y<ZB{A) YCB(A) 

where N^{Y) = |A| - \S{X) U dA\ + C{X\dA). 



2.2 Low temperature expansion of the bulk partition 
function 

We give in this subsection an expansion of the partition function Q°{A) at 
"temperature" /3b > ^. The expansion is mainly based on a duality property 
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and we first recall geometrical results on Poincare and Alexander duality (see 

e.g. PB],[Il,[I3],[in|). 

We first consider the lattice Zf- and the associated cell-complex L whose 
objects Sp are called p-cells (0 < p < d): 0-cells are vertices, 1-cells 
are bonds, 2-cells are plaquettes etc.: a p-cell may be represented as 
(x; CTiCi, CTpCp) where x G Z"^, (ci, e^) is an orthonormal base of and 
a a = ±1, a = 1, ...,d. Consider also the dual lattice 

(Z'^)* = |x = (xi + ^, Xrf + ^) : Xc, e Z, a = 1, 

and the associated cell complex L*. There is a one to-one correspondence 

Sp ^ 4^p (2.5) 

between p-cells of the complex L and the d — p-cells of L*. In particular to 
each bond si corresponds the hypercube s^_^ that crosses si in its middle. 
The dual E* of a subset E' C L is the subset of element of L* that are in the 
one-to-one correspondence ()2.5|) with the elements of E. 

We now turn to the Alexander duality in the particular case under con- 
sideration in this paper. Let Y C B{A) be a set of bonds. We define the 
A-dual of Y as 

Y = {B{A)\Yy (2.6) 
As a property of Alexander duality one has 



Y 



\B{A)\Y\ (2.7) 



ArX(F) = N^iY) (2.8) 

where Nci{Y) denote the number of independent closed {d — l)-surfaces of 
Y. We thus get 



g°(A) = g^'-l^^^)! Yl q-^"^^^^""^'^^^ (2.9) 

Yc[B(A)r 



This system can be described by a gas of polymers interacting through 
hard core exclusion potential. Indeed, we introduce polymers as connected 
subsets (in the R"' sense) of {d — l)-cells of L* and let V{A) denote the set of 
polymers whose {d — l)-cells belong to [5(A)]*. Two polymers 71 and 72 are 
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compatible (we will write 71 ~ 72) if they do not intersect and incompatible 
otherwise (we will write 71 72). A family of polymers is said compatible 
if any two polymers of the family are compatible and we will use P(A) to 
denote the set of compatible families of polymers 7 G 'P(A). Introducing the 
activity of polymers by 

<^„(7) = q-Mi\+NAi) (2.10) 

one has: 

g°(A) = ^ n^oW (2.11) 

yeP(A) 7er 

with the sum running over compatible families of polymers including the 
empty-set with weight equal to 1. 

We will now introduce multi-indexes in order to write the logarithm of this 
partition function as a sum over these multi-indexes (see j2HI)- A multi-index 
C is a function from the set 'P(A) into the set of non negative integers, and 
we let suppC = {76 'P(A) : C(7) > 1}. We define the truncated functional 

*.(c) = nS^n^oW"'^' (2-12) 

where the factor a{C) is a combinatoric factor defined in terms of the connec- 
tivity properties of the graph G{C) with vertices corresponding to 7 G supp C 
(there are C{^) vertices for each 7 G suppC ) that are connected by an edge 
whenever the corresponding polymers are incompatible). Namely, a{C) = 
and hence $o(C') = unless G{C) is a connected graph in which case C is 
called a cluster and 

a{C)= J2 (-1)''^^^' (2.13) 

GCG{C) 

Here the sum goes over connected subgraphs G whose vertices coincide with 
the vertices of G{C) and |e(G')| is the number of edges of the graph G. If the 
cluster C contains only one polymer, then 0(7) = 1. In other words, the set 
of all cells of polymers belonging to a cluster C is connected. The support of 
a cluster is thus a polymer and it is then convenient to define the following 
new truncated functional 

$(7)= J2 '^o(C) (2.14) 

Cisupp C=-f 

As proved in paper I, we have the following 
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Theorem 2.1 Assume that > 1/d and Cou^q < 1^ where Vd 

^224(d-i)^ and Co = [l + 2'^~'^{l + VlT^)] exp 



QO(A) = g/3,|B(A)| J I (2.15) 

with a sum running over (non-empty) polymers, and the truncated functional 
$ satisfies the estimates 

|$(7)l < bl fcoz/dg-^''+^)'''' (2.16) 



The proof uses that the activities satisfy the bound (^0(7) < q-iPb-^/d)\'y\ 
(because iVci(7) < \l\ld) and the standard cluster expansion. The details 
are given in Ref. |llj . 

2.3 Hydra model 

We now turn to the partition function Z^iVL). We will, as in the previous sub- 
section, apply Alexander duality. It will turn out that the ratio ZP(fi)/Q°(A) 
of the partition functions entering in the definition ()1.4|) of the surface free 
energy can be expressed as a partition function of geometrical objects to 
be called hydras. 

Namely, we define the A-dual of a set of bonds X C B{VL) as 

x = {B{n)\xy (2.17) 

This transformation can be analogously defined in terms of the occupation 
numbers 

^ ^ f 1 if & e X ^2 18) 

1^ otherwise 

For a configuration n = {nb}beB{n) C {0, l}^^^^ we associate the configura- 
tions n = {nsjseiBm* C {0, 1)1^(^)1* given by 

Hfe. = 1-^5, beB{n) (2.19) 

where b* is the {d — l)-cell dual of b under the correspondence ()2.5|) : (see 
Figure 2). 
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FIGURE 2: A configuration X (full lines) and its A-dual X (dashed lines). 



For any set of cells X we will use the decomposition X = U ZhU Yf, 
where Xg is the set of cells whose dual are bonds with two endpoints on the 
boundary surface E, Zh is the set of cells whose dual are bonds with one 
endpoint on the boundary surface S and one endpoint in the bulk A and the 
remaining Yfy is the set of cells whose dual are bonds with two endpoints in 
the bulk A. Thus, for the decomposition X = XsU Xh introduced above, we 
have 



+ 



Xs 



\B{Q)\B{^)\-\X,\ 



We let Bq be the set of bonds that have an endpoint on the boundary layer 
Lq and the other endpoint on the layer L_i = {i E IL \ = —1} and let 
Nci{X) be the number of independent closed surfaces of X U B^: Nci{X) = 
Nci{X U Bq). As a result of Alexander duality, one has 

Ar°(X) = iVei(X) 

Denoting by Bi{Q) the set bonds that have an endpoint in dg^l the other 
endpoint in L \ fi, we have furthermore 

N^\X) = N,,{XU[B^{Q)]*) 

These formula lead to the following expression for the partition function ()2.3p 



zp{n) = q 



A|B(S)|+/36|B(n)\B(E)| 



-l3s\Xs\-Pi{\zA + \Yf,\)+N,i{X)+xl 
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where = and Xn = N^ii^ U - A^ci(^)- Notice that this 

Boltzmann weight equals the Bohzmann weight of the bulk partition function 
(im for those X C [B{A)]* i.e. if X, = and Zi, = 0. They can thus be 
factorized in the ratio of the two partition functions. Namely, we define 
hydras as components of {d — l)-cells not completely included in [i?(A)]*. 

Definition 2.2 A connected set of (d—l) -cells S C (intheW^ sense) 

is called hydra in Vt, if it contains a cell whose dual is a bond with at least 
one endpoint on the boundary surface S. 

Definition 2.3 Given an hydra 6 C [-B(fi)]*, the components of 6 included 
in [B(I1)]* are called legs of the hydra, the components included in [B{A)]* 
are called heads of the hydra and the remaining components are called bodies 
of the hydra. 

The dual cells of bodies of hydras are bonds between the boundary layer and 
the first layer hi = {i & h \ id = 1}; (see Figure 3). 



FIGURE 3: A hydra, in two dimensions (a dimension not considered in this paper), 
with 5 feet (components of fuU lines), 2 bodies (components of dashed hnes), and 3 heads 

(components of dotted hnes). 

We let TC{fl) denote the set of hydras in Q. Two hydras 6i and ^2 are 
said compatible ( we will write 6i ^ 62 ) if they do not intersect. A family 
of hydras is said compatible if any two hydras of the family are compatible 
and we let H(f2) denote the set of compatible families of hydras 6 G H{fl). 

Clearly, a connected subset of cells included in [i?(f2)]* is either a hydra 
6 G T-C{Q) or a polymer 7 G 'P(A) (defined in Subsection 12. 2^ . Therefore 
any subset of [-B(f2)]* is a disjoint union of a compatible family of hydras 
X G H(fi) with a compatible family of polymers Y G P(A). 
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The partition function Zp{Q) given by (j2.H) reads thus: 
ZP^n) = qMB{i^)\+MB{n)\B{j:)\ ^2.20) 
X q-MX4-t3U\Zb\ + m)+N,i(x)+x^, ^ JJ 

xeu{n) y6P(A):y~x 7ey 

where the compatibihty Y ^ X means no component of Y is connected with 
a component of X. 

According to Subsection 123 the last sum in the RHS of the above formula 
can be exponentiated as: exp ^J^jeviAy-y^x ^i^/)^- Hence dividing the above 
partition function by the partition function Q°{A) we get by taking into 
account Theorem 12.11 

zp{n) 



SP(Jl) = ^^^^^ = ^A|B{S)|+/36(|B(f^)h|B(S)h|B(A)|) (^2.21) 



X q-MXs\-M\Zt\ + \Y,\)+N,i{X)+xf,^^^ ) ^ ^^^^ 

Hereafter the incompatibility ^ X means that no component of X is 
connected with 5. 

HP(fi) is thus the partition function of a gas of hydras X = {6i, . . . , 5„} 
interacting through hard-core exclusion potential and through a long range 
interaction potential (decaying exponentially in the distance under the hy- 
pothesis of Theorem I2.ip defined on the polymers of the bulk. 

If we neglect this long range potential, and if we moreover restrict to 
consider only hydras without head, the system of hydras will reduce itself 
to a (c? — 1) Potts model with two-body interaction coupling K and mag- 
netic field J (i.e. with formal Hamiltonian H = —^(^t j'^d^^ K6{cri,aj) — 
k) ieLo fceLi J^i^ii 0) )■ This model undergoes a temperature driven first 
order phase transition, whenever q is large enough and d > 3 [Sj. We will 
show that it is also the case for the hydra model 1)2.211) implementing the fact 
that the heads of hydras modify only weakly their activities and that the long 
range interaction potential decays exponentially (the needed assumptions are 
close to those of Theorem 12.1)) . To this end it is convenient to first rewrite 
this potential in terms of a model of aggregates. 

Let us introduce the (real-valued) functional 

^(7) = e-*(^) - 1 (2.22) 
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defined on polymers 7 G -P(A). An aggregate A is a family of polymers 
whose support, suppA = U^g^7, is connected. Two aggregates Ai and A2 
are said compatible if supp Ai fl supp = 0. A family of aggregates is said 
compatible if any two aggregates of the family are compatible and we will 
use A (A) to denote the set of compatible families of aggregates. Introducing 
the statistical weight of aggregates by 



uiA) = n *(7) (2.23) 

■y&A 



we then get: 



exp ^ - J2 '^(7) \= n (1 + "^(7)) 

= E n ^(^) (2.24) 

^eA{A) A€A;AooX 

where A X means that every polymer of the aggregate A is incompatible 
with X. Since the support of aggregates is a connected set of {d — l)-cells, 
i.e. a polymer, it is convenient (as it was done for clusters in Subsection 2.3) 
to sum the statistical weights ()2.23|) over aggregates with same support. We 
thus define the weight 

^(7)= E ^(^) (2.25) 

j4:supp y4=7 

with A X, to get 

;=;P(^) ^ ^A|B(E)|+/36(|i3{C)|-|iJ(S)|-|B(A)|) ^2.26) 

X q-MXs\-M\Zi\+m)+N,i{x)+xl ^ J-j- ^^^-^ 

The system is thus described by two families: a compatible family of 
hydras (a subset of [i?(L)]*) and a compatible family of polymers (a subset 
of [i?(L\Lo)]*) each of these polymers being incompatible with the family 
of hydras. 

We will define in the next subsection the diluted partition functions for 
our system. This partition function differs only from the "physical" partition 
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function (|2.26|) by a boundary term and thus both partition functions lead to 
the same free energy. The recurrence relations of Lemma f3. II below, allow to 
expand the diluted partition functions in term of matching signed contours 
interacting through hard-core exclusion potential. 

2.4 Diluted partition functions 

Note first that even though our model is defined on a (i-dimensional box Q 
it has a (d — l)-dimensional structure and the highest order of the logarithm 
of partition functions behaves like 0(|S|). It will be convenient to consider 
f2 as a set of lines and its dual Q* as a set of columns. 

We let a line L{x) be a cylinder set of sites of L whose projection on the 
boundary layer is the site x and whose height is less than a given number 
M: L{x) = {i eh [ii, id-i) = x eho, id < M}. We let Lm be the set of 
all such lines. The dual of a line is called column and is thus a set of (i-cells 
of the complex L*. For Q C Ljv/, we let S = ^7 fl Lg, be its projection on the 
boundary layer, A = 1] \ S and = |S| be the number of columns of Q* 
(or of lines of Q). 

Consider a site x G L and its dual d-cell x*. We shall use S{x*) to 
denote the set of {d — l)-cells of the boundary of x* (there are the dual 
cells of the bonds whose x is an endpoint). For a set of (i-cells D, we let 
S{D) = Ux*(=dS{x*) be the union of the boundaries of the d-cells of D. 

Next, it can easily be checked that the configuration {X° = 0, F = 0) 
and the configuration (X°^ = [5(Lo)\-B(L\Lq)]* , F = 0) are ground states 
of the system. 

We will use H^IQ) to denote the set of compatible families of hydras 
defined on£{n*)n [B(L)]* that coincide with Xp on£{[dn]*), and use P^'^ (A) 
to denote the compatible families of polymers defined on S{Q*)\{S{T,*) U 

s{[dn]*). 

For such configurations the Boltzmann weight in ()2.26|) reads 

^-MXs\-M\Zt\+\n\)+N,i{x) 

since for those X G H°^(fi) one has N^i{X U [Bi{n)]*) = N^i{X). 

We define, for (any) volume Q Gh, the diluted Hamiltonian of a config- 
uration X = XP 

H^\X) = e.*(X) - iVci(X n £{n*)) (2.27) 
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where the energy per cell is defined by 



2 



for the (i-cells of the surface and by 



Zbn£{x*) 



if X G Ln 



if X G LXLn 



for the (i-cells of the bulk. 

The diluted partition function is defined by 



(2.28) 



Up to a boundary term 0(9S) one has InS^ (Q) = [{d — l)Ps + Pb] ||^|| lng + 
InEf^ (n), hence 



hm -— In Ef' (n) =g,+ [{d -l)f3, + (3t,] In q 



(2.29) 



where f2 | L means that we take first the limit M — > oo and then the limit 
S t Lq in the van- Hove or Fisher sense [SI]. 

Notice that the diluted Hamiltonian on ground states reads on set of 
columns Q* C L^^: 

H^\XP) = ep\\n*\\ (2.30) 



with the specific energies 



Co = 

Cof = {d-l)p, + pb-l 



(2.31) 



3 Surface transition in the bulk low temper- 
ature regime 

3.1 Contours and Peierls estimates 

We first define the contours of our system. 

Let fl C Lm, ^* its dual set and (X, Y) be a configuration of our system 

inn-. X e Hp(^]), Y G P'^^{n), f ^ x. 
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A d-ceW X* G Q* is called p-correct, if X coincides with the ground state 
on the {d — l)-cells of the boundary S{x*) of x* and the intersection 

Y n S{x*) = A column is called p-correct if all the d-cells of the column 

are p-correct. 

Columns and d-cells that are not p-correct are called incorrect. 

The set of incorrect columns of a configuration {X, Y) is called boundary 
of the configuration {X, Y). 

A triplet F = {supp F, X{r), Y{r)}, where supp F is a maximal connected 
subset of the boundary of the configuration (X, Y) called support of F, X(F) 
the restriction of X to the boundary £^(suppF) of the support of F, and 
Y{T) the restriction of Y to £(supp F), is called contour of the configuration 
{X,Y). Hereafter a set of d-cells is called connected if the graph that joins 
all the dual sites i,j of this set with d{i,j) < 1 is connected. 

A triplet F = {supp F, X(F), F(F)}, where suppF is a connected set of 
columns is called contour if there exists a configuration (X, Y) such that F 
is a contour of {X,Y). We will use |F| to denote the number of incorrect 
cells of suppF and ||F|| to denote the number of columns of suppF. 

Consider the configuration having F as unique contour; it will be denoted 
{X^,Y^). Let Lp{r) be the set of p-correct columns of \ suppF. Ob- 
viously, either a component of Lo{T) is infinite or a component of I/of(F) is 
infinite. In the first case F is called contour of the ordered class or o-contour 
and in the second case it is called of-contour. When F is a p-contour (we 
will let F^ denote such contours) we use Ext F to denote the unique infinite 
component of Lp(r); this component is called exterior of the contour. The 
set of remaining components of Lp{T) is denoted IntpF and the set Lm^piV) 
is denoted Int^F. The union IntF = IntfF U IntfoF is called interior of the 
contour and V{T) = supp F U IntF. 

Two contours Fi and F2 are said compatible if the union of their supports 
is not connected. They are mutually compatible external contours if V^(Fi) C 
ExtFs and V{T2) C ExtFi. 

We will use G'(F^) to denote the set of configurations having F^ as unique 
external contour. The crystal partition function is then defined by : 

E ?"''"^^'^^^n^(^) (3-1) 

{X,Y)eG(TP) 765^ 
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Lemma 3.1 The following set of recurrence equations holds: 



-ep\\Ext\\ 



(3.2) 



i=l 



Here the sum is over families {F^, . . . , T^}ext of mutually compatible external 
contours in Q (suppT^ C flmt = {i & fl : d{i, La/ \ fi) > 1} ), \\Ext\\ = — 
^||V(r^)|| where ||V(r^)|| is the number of columns o/V"(r^); 



m£{o,o/} 



where: 



supp L f ^ ' 



n 



Proof. We have only to observe that for any X G 11^(0) 



(3.3) 



(3.4) 



(3.5) 



where the sum is over all contours of the boundary of the configuration 
(X, y = 0) and Ly{}C) fl Vl* is the number of p-correct columns inside 
of this configuration. ■ 

Lemma f3. II gives the following expansion for the partition function 



L-^ ]^ ! ■ ■ ■ j-^ n -T comp 



(3.6) 



where the sum is now over families of compatibles contours of the same class 
and 

|S5^i(Int„r^') 



^(rf) = Q{Y^)q'^\\^''\\ 



sf(int„rp) 



(3.7) 



where ||r^|| is the number is the number of columns of supp F*' and m ^ p. 

To control the behavior of our system, we need to show Peierls condition, 
that means that Q{V'P)q^p\\^^\'^ has good decaying properties with respect to 
the number of incorrect cells of supp F^*. We use in fact the modified Peierls 
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condition introduced in Ref. [21] where g{TP)q^p^^^''^^ is replaced by g{TP)q-^^^^^^ 
with e = min (cq, Cof). 
Let 



g-r _ /'2(3d-2)^-2i^ ^ 3c2"'+V|g^-^'') (3.8) 



1 - Qcu^q 

where c = 8e(e — 1)cq and = ci^2'*('^~^\ We have the following 



Proposition 3.2 Let S C L^,^ be a finite connected set of columns, assume 
that 2 < Pb < ^ and 6cz/|g~3+/^b <- ^/^g^ fg^. q// ^ 

J2 |^(mg-"^'" <e-^i'^l' (3.9) 
rp:sMpprp=s 

where \\S\\ is the number of columns of S. 
The proof is postponed to the Appendix. 

The recurrence equations of Lemma 13.11 together with the Peierls esti- 
mates ()3.9p allow to study the states invariant under horizontal translation 
(HTIS) of the hydra system as in paper I. This is the subject of the next 
subsection. 

3.2 Diagram of horizontal translation invariant states 

To state our result, we first define the functional 

K,{S)= (3.10) 

rPisupp TP=S 

Consider the partition function Sp'^(r2) ()3.6|) and for a compatible family 
{V\, ^n}co'[Dp p-contours, denote by 5*1, ...,Sn their respective supports. 
By summing over all contours with the same support this partition function 
can be written as the partition function of a gas of polymers S with activity 

Kp{S) = J2 ^(X^) interacting through hard-core exclusion potential: 

rp:supp rp=s 

n 

Ef\n) = q~^^\l^l\ ^ n^^(^^) (3-11) 
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Here {Si, •S'^j^.^^p denotes compatible families of polymers, that is 
d{S*, S*) > 1 for every two Si and Sj in the family: recall that by definitions 
of contours a polymer 5* is a set of columns whose graph that joins all the 
points of the dual of the columns of S at distance d{i,j) < 1 is connected. 

Next, we introduce the so-called truncated contour models defined with 
the help of the following 

Definition 3.3 A truncated contour functional is defined as 

P [0 otherwise 

where \\Kp{S)\\ = J2rp-supprp=s ^'^^ a > is some positive parameter 

to be chosen later (see Theorem \3.5\ below). 

Definition 3.4 The collection {S,p} of all p- contours T^with support 
suppT^ = S is called stable if 

<e-"ll^il (3.13) 

z.e. zfK,iS)=K;iS). 

We define the truncated partition function Sp(f2) as the partition function 
obtained from Sp''(f2) by leaving out unstable collections of contours, namely 

n 

E'^(Q)=q-^.m Yl^m (3.14) 

comp 

n 

= ^ 1[k;{S.) (3.15) 

Here the sum goes over compatible families of stable collections of con- 
tours. Let 

be the metastable free energy of the truncated partition function 

For a large enough, the thermodynamic limit (j3.16|) can be controlled by 
a convergent cluster expansion. We conclude the existence of hp, together 
with the bounds 

g-Ke-^ia^ni ^ E' {n)e^^^^^^^ < (.'^'^'"I'^M (3-17) 
\hp — Cpliaql < Ke"" (3.18) 
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where k = HdWY where Kd = ^+^ e^+i is the cluster constant [20] and 
= ?)^~^ — 1; = dVt n Lq in the way defined in Subsection 11.21 



Theorem 3.5 Assume that 1/d < Pb < 1 and q is large enough so that 
e-" = e- 
such that 



-t+2k'+3 ^ ^^gj^ there exists a unique jSl = ^(1 — Pb) + 0{e 



(i) forPs = Pl 



for both boundary conditions p =o and p =of, and the free energy of 
the hydra model is given by gf+ [{d — l)/5s + Pb] In q = ho = hof 

(ii) for (3s > (31 

and go + [{d - l)(3s + (3b\ Inq = ho < hof 

(iii) for (3s < (31 

and gof+ [{d - l)Ps + Pb] Ing = hfo < hj 

Proof. Starting from the Peierls estimates given in Proposition 13.21 and 
the definitions of this subsection the proof is the same as that of Theorem 3.5 
in paper I. Let us only recall that to exponentiate the partition function 
Zp{iT) = g^pll^llSp(fi) we define the truncated functional associated to K'^ 

^^ix) = YfYlsyUKis)''''' (3.19) 

defined on the multi-indexes X associated to the polymers (a multi-index 
being a function from the set of polymers into the set of non negative integers, 
and a{X) is defined as in ()2.13|1 ). The number of polymers S with number 
of columns US'!! = n and containing a given column can be bounded by z/" 
where u = (3^^^^ — 1)^ as in paper I: this is because the chosen definition for 
connectedness of columns here is the same as that for connectedness of lines 
in paper I. 

As a result of the standard cluster expansion [2011211, we get for kc""' < 1 



Z,(fi)=exp|5^<|.^(X)| 
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where the sum is over multi- indexes whose support suppX ={S':X(S')>1} 
belongs to fl. The series J^x snvvx^L |'^"^(^)| absolutely convergent and 



^Jf isuppX 

satisfies the bound 

X.snppX 3L 



<l>^(X)|<Ke"" (3.20) 



Let us introduce the Gibbs states ( ■ associated to the boundary con- 
ditions p G {f;of}. Theorem above show also that at Ps = Pi 

(n,)° = 1-O(e-0 (3.21) 
(nfe)°f = O(e-0 (3.22) 

for any bond b of the boundary layer and any bond b between the boundary 
layer and the first layer. 

Indeed by the correspondence ()2.19|) these equations are equivalent to 

= 0{e~^) (3.23) 
(n,.)°^= 1-O(e-0 (3.24) 

for the dual cells of the bonds under consideration. By definition of contours, 
with ordered (o) boundary conditions, any such cells are occupied only if 
there is an ordered contour surrounding it and that the correlation functions 
are controlled by the contour model cluster expansion. With ordered-free 
(of) boundary conditions, such cells are empty only if there is a of-contour 
surrounding it and again the correlations are controlled by cluster expansion. 
Obviously, the relation (j3.23|) holds true for any Ps > Pi while the relation 
fHT^ hold true for any Ps < Pi- 

This shows in particular that the derivative -^Qo of the free energy Qo 
with respect to the surface coupling constant K is discontinuous near K = 
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Note added after publication: Let us mention Refs. P in which, besides 
numerical results analytical calculations in the large-q limit, mainly in 2d, 
are presented. 



Appendix: Proof of Proposition 13.2 



We begin the proof by considering contours F = |supp F, X'", y^^j (where 
jx'", ^'"l is the configuration having F as unique contour) without polymers, 

i.e. = We have the decomposition X^ = Xf U U where Xf = 
Xr n [5(Lo)]*, Y^ = X^n [5(L \ h,)Y, and Z[ = \ (Xf U Y,^). 

A d-cell X* G Lq will be called regular if the {d — l)-cells of its boundary 
that belong to the boundary layer are either all empty or all occupied. It 
will be called irregular otherwise. We denote by -Ro(r) the set of correct d- 
cells of the contour F: RoiT) = {x* e Lq : jX^ n^(a;*)| equals or 2rf - 1 } . 
We let /o(r) be the set of incorrect rf-cells of the contour F : /o(F) = 
{x* e Lo : 1 < \X^n£{x*)\ < 2(rf- 1)}. 



Lemma A.l 



^(Dg^liril < q-^\Hn-{P^-m\ (A.l) 
Proof. By Lemma (3. II and definition ()2.27|) . one has 

g(jp-^ ^ qT,^''eRo(r)uio{r)''^''i^)-(^b\Yt\+N^iiX) 



where e^*(X) = -f X, n £{x*) ~ (3b Z^n £{x*) 



and to simplify formulae 



we put hereafter X, X^, and Y^ instead of X^, Xf , Z^ and Y^ . We define 

iVei(x, u Zfc) = iVei(x, u z,) - N,,{xns{x*)) 

as the number of independent closed surfaces that are not boundaries of an 
occupied d-cell of the surface. This leads to the decomposition 

f,(r)g^liril = g-A4r)-B.(r)-A,(r^) (^.2) 
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where 



As{T)= [e.*iX)-e-N,i{Xn£ix*)) 

a;*ei?o(r) 

B,(r)= [e:.'W-e 
a;*e/o(r) 

Afe(r) = 13b Inl - fiVd(X) - iVd(X, U Z,) 



- iVel(X, U Zb 



Clearly 



(A.3) 
(A.4) 
(A.5) 

(A.6) 



A,(r) > 

Indeed the regular d-cells are either empty in which case the term inside 
brackets in ()A.3|) equals 



-e = - min {cq = 0; Cof = (d- l)/^^ + - 1} 

or they are occupied in which case this term equals 

id - + /^fe - 1 - min {co = 0; Cof ={d- + /Jf, - 1} 

Let us now bound Bs(r). We first notice that for incorrect (i-cells x* of 
the surface, 



ex*(A) - e = — 



> 



Xs n £{x* 



+ /3b 



ZhHSix* 



1-A 
rf- 1 



+ (3b 



Zh n ^Cx* 



2(rf- 1) 

Furthermore the number Nc\{Xs U Zh) may be bounded as 



(A.7) 



Ael(A,UZfc) < 



E*:ii-l<|Xn£(a:*)|<2{d-l)-l 



2\xn£{x'')\ 

22('i-l) 



< 



E 



a:*:ti-l<|if n£(2:*)|<2{d-l)~l 
|Zjn£(a;*)|=l 



1 

2 



If an incorrect site of the surface is such that 



z,, n Six*] 



(A.8) 
0, then 



necessarily 



X, n Six* 



1 and thus such site gives a contribution 

2(rf- 1) 
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toBs(r). Let us now consider those incorrect (i-cells for which Zbr\S{x*) 
1. Starting from ()A.7|) we get for such cells 



e...(X) - e > (3b + 



2{d- 1 
1-/36 



n Six*] 



> 



2{d-l) 
2{d-l) d 



2{d- 


1) 


+ 1 - 


Xsns{x*) 




x^nsix*) 


- 1 






2{d- 


1) 


- + 


2{d-l 


) 


2{d- 


1) 


+ 1 - 


X, nf(x*) 


+ - 


Xsn£{x*) 


- 1 



Pb 



x.^nsix*] 



+ 1 



2{d-l) 2d 
where for the second inequality we take into account that 



x, n Six*] 



< 



2{d — 1) and I3b>l/d. When furthermore 



x.r\£ix*] 



> d — 1, one infer 



I -13b 



2(d-l) 2 

and it thus follows from ()A.4jl . ()A.8jl . and ()A.7j) . that each incorrect cell with 



Z^nEix* 



1, gives also a contribution at least 2(d-i) ^^(r). Therefore 



<q 



■2^l^o(r)| 



Consider finally, the quantity Ab(r). We will prove the inequality 

Yb 



d 



(A.9) 



(A. 10) 



X,i(X)-X,i(X,UZfe) < 

Notice first the obvious inequality 

iVei(X) - iVd(X, u z,) < Xd(n U 5*J 

where Bqi is the set of bonds between the boundary layer and the first layer 
and Nc\iY}y U B^^) is the number of independent closed surfaces of Yf, U -Bq^. 



The number 



can be written 



Y. 



G [LXLof : I YbD^Cx* ) I > 1 



Y, n £ix* 
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(because each {d — l)-cell belongs to the boundary of two (i-cells). 

For the configurations Y/, that do not hit Bq^ (meaning that there are no 
cell of Yb connected with Bq^ in the M.'^ sense) we get immediately ()A.10|) as 
already used in the proof of Theorem 12.11 closed surfaces of minimal area 
are d-cells and the number of {d — l)-cells in the boundary of d-ceW equals 
2d. 

For the configurations Y^ that do hit -Bq^, notice first that can be 
written 



Yh 



E 

x*e[L\Lo]*:|Yi,nf (a;')|>2 



Yh n Six* 



E 

x*e[L\Lor:|Y6n£'(a;*)|=l 



Y,ns(x*] 



Since Yh is finite, the set of cells with 
furthermore for any [d — l)-cell of -Bgi t 



Y, n Six*] 



(A.ll) 

1 is non empty and 



lere exists a d-cell x* above it such 



that 



YhHSix* 



= 1 that can not contribute to Nci{Yi, U -Bqi). Now closed 

surfaces of minimal d-cells having in their boundary 2d — 1 cells of 

Yh and a cell of B^^. Thus for such surfaces we have a contribution d — 1/2 
coming from the first term of the RHS of ()A.11|) and a contribution 1/2 
coming from the second term. 
This implies flA.lOj) giving 



-A6(r) 



<q~ 



(A.12) 



which in turn implies the lemma by taking into account ()A.2|) . ()A.6|) and 

Considering still contours F = jsuppF, X^, F'"} without polymers, i.e. 
{Y^ = 0) we have the 

Lemma A. 2 Assume that jd^ > h, and 2i'dq^~^^ < 1, then 



f^(F)g^liril < [2(''~'k'^- + 2 

r:suppr=s 



d+l 



\\s\\ 



1 



2udq^-f^^ 

(A.13) 



which shows that, whenever q is large enough, the Peierls condition holds 
true for the class of contours without polymers. 
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Proof. 

First, observe that for contours F with support supp T = S and number 



of irregular cells of the boundary layer |/o(r)| = k one has Yi, 
\Sm\ > \\S\\-k. Therefore, 

r:suppr=5 o<fe<||S|i rp:|/o(r)|=fc 



< J2 /'ll'5|l'\2(2<i-i)fe2ll^ll-V^' (A.14) 

n^u^waw V ^ / 



0<fe<|| „ 

m 

n<2||5|| Si3''l,-;Sn3sn si,...,s„ 

l«ll + ... + l'5nl>llS||-fc SaeS;SQ*S*j^ 



X 



Here the binomial coefficient ("f") bounds the choice of irregular cells of 
the dual of the boundary layer while the factor 2(^'^^^)'^'2'l'^ll^'^ bounds the 
numbers of contours with H^*!! columns and k irregular cells; the notation 
Sa -Sqi means that a, {d — 2)-cell of the boundary of the {d — l)-cell 
belongs to the boundary 8(3^-^). Then 



E ^(F)g^liril < Yl ('l'^ll]2(2'^-^)^2ll^ll-^g-^'= (A.15) 

r:suppr=5 0<A:<||S|| ^ ^ 

X E C^-:'"^") E nM-)"' 

n<2||5|l ^ ^ mi+...+m„>||5||-A:i=l 

Here the binomial coefficient ('''^ 1^^""^") bounds the choice for the components 
Si, ...,5n of Yf, to hit the boundary layer at Si, ...,s„. The above inequality 
yields 



i)k2\\s\\-k -^^k 



r:suppr=5 0<fe<||S|| ^ ^ 

E E (-^«-^')" 

n<2||S|| ^ ^ m>||5||-A: 



:supp 

X 



0<ifc<||5|| \ ^ / 



X (2udq-y-^-f''')) ^ (A.16) 

^ ^ 1 — 2i/^0d 
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that gives the inequahty of the lemma. ■ 

We now turn to the general case of contours with non empty polymers 
and first give a bound on the activity (7) of polymers. 

Lemma A. 3 Assume that Pb > and cv'j^q~^~^^ < 1 with c = 8e(e — l)co 
and Vd = (2(i)^, then 

1^(7)1 < (ci^dg^"^')'^' (A. 17) 

Proof. Let us first recall the definition ()2.25|) : ^'(t) = I^AsuppA=7'^(^) 
where the weights of aggregates are defined by (see ()2.22|) and ()2.23|1 ): uj{A) = 

n^g^(e-*(^)-l). By TheoremOwe know that |<l>(7)| < (^ecoi^dg^"^") (< 
1) for q large enough. Since for any |a;| < 1, |e^^ — 1| < (e — 1) we have 

1^(7)1 = |e-*(^)-l| < (e-l)|$(7)| < ((e-l)ecoz/dg^-^'')' 



l7l , , 



(A.18) 



Then, 



A:suppA=7 n>l 71,...,7n: j=l 

supp {71.---.7n}=7 



n 



< 5^ 21-^1 Yi n^""'"''' 

n>l 7l331>--->7n3sn: j = l 

supp {7i,...,7n}=7 

n>l mi,...,m„: j=i 

mi + ...+m,i>|7| 

n 

„>1 mi,...,m„: j=i 
mi+...+m„>|7| 

Here, we used as in the proof of Theorem 12.11 that the number of polymers 
of length m containing a given bond or a given vertex is less than z/^; the 
term 2^^^^ bounds the combinatoric choice of the cells Sj G 7j, because 7 being 
connected, it contains n — 1 such intersecting cells (see [T3]). 
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We put k = nil + ... + and notice that there are at most such 
numbers to get 

A:suppA=7 l<?i<A:fc>|7| ^ ^ 



<5:(4..e-)'=EG 



2 



fc>|7| 



fc>|7| 



< 



l7l 



(A.20) 



provided that ^i^jl'' < 1- The lemma then follows by assuming that 

We finally turn to the 
Proof of Proposition 13.21 

Consider a contour T = |supp F, X'", F'"! and as above the decomposi- 
tion = X[ U Zl U Y^. Consider also the union fb = U Y^. Notice 
that the set T = U U Tj, is a family of hydras and there are at most 
sl"^*! contours corresponding to this family: this is because a {d — l)-cell in 



Th may be occupied either by Yj^ or by Y^ or by both. Let 



(A.21) 



The above remark on the number of contours associated to T and Lemma FA .31 
imply 



eliril 



|J"o(r)| 



< q 3sup <^ qd 



IMEli 



(A.22) 



The rest of the proof is then analog to that of Lemma IA.2I starting from 
Lemma IA. 31 and replacing qd-l^^ by 3cujqd~'^>' . It gives 

J2 l^(r)| g^ll^ll < (2^^'-%-^ + 3c2^+VJg^-/^^)'"" 

(A.23) 

provided 6cz/Jgd~^^ < 1 and ends the proof of the proposition. ■ 
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